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We have used small-angle-neutron-scattering (SANS) and ac magnetic susceptibility to investigate 
the global magnetic field H vs temperature T phase diagram of a single crystal Nb in which a first- 
order transition of Bragg-glass melting (disordering), a peak effect, and surface superconductivity 
are all observable. It was found that the disappearance of the peak effect is directly related to a 
multicritical behavior in the Bragg-glass transition. Four characteristic phase boundary lines have 
been identified on the H — T plane: a first-order line at high fields, a mean-field-like continuous 
transition line at low fields, and two continuous transition line associated with the onset of surface 
and bulk superconductivity. All four lines are found to meet at a multicritical point. 



PACS numbers: 74.25.Qt 61.12.Ex 

An outstanding question concerning the Abrikosov 
vortex state of type-II superconductors is whether a gen- 
uine order-disorder transition can still occur in vortex 
matter even though true crystalline order cannot be at- 
tained due to random pinning by impurities Jl| . There 
are convincing theoretical arguments 0, HI 0, U and nu- 
merical evidence |(| suggesting that, instead of a true 
vortex crystal, a novel Bragg glass phase with quasi-long- 
range order can exist in bulk samples with weak random 
pinning; hence a true order-disorder transition can oc- 
cur when the topological order of the Bragg glass is de- 
stroyed, by thermal fluctuations and/or random pinning. 

However, it is still controversial as to whether a Bragg 
glass melting (or disordering) transition is the underlying 
mechanism of the well-known anomaly of "peak effect" 
in weak-pinning type-II superconductors • Recent neu- 
tron scattering experiments on Nb pi , and VsSi as 
well as STM studies of 2H-NbSe 2 H3 all suggested a 
disordering transition at the peak effect, and the phase 
transition appears to be first order p|. However, it has 
been reported that some samples of similar quality, e.g., 
having only weak pointlike pinning centers, do not show a 
peak effect, nor a disordering phase transition ^l|. This 
raises an obvious, but intriguing, question: Is the fate of 
the peak effect, i.e. appearing or disappearing, related 
to a multicritical behavior in the phase transition into 
the Bragg glass? In this Letter, we report the first di- 
rect evidence that the disappearance of the peak effect is 
related to a multicritical point on the Bragg-glass phase 
boundary. 

Our experiment was carried out using the 30m SANS 
instruments NG7 and NG3 at the NIST Center for Neu- 
tron Research on a Nb single crystal (99.998% in purity) 
in which both the peak effect and the first-order Bragg- 
glass melting (disordering) transition were observed at 



the same temperatures p|- The sample has a zero-field 
T c — 9.16 K, and an estimated Ginzburg-Landau param- 
eter Ki(0) = 2.0. The mean wavelength of the incident 
neutron beam was A = 6.0 A and the wavelength spread 
11 % (FWHM). The experimental configuration is shown 
in the inset of Fig.l (a). A cadmium mask was used such 
that only the central portion of the sample was exposed 
to the incoming neutron beam. The scattered neutrons 
were captured by a 2D detector of 128x128 pixels (the 
pixel size is 0.5 cm by 0.5 cm) 15.3 m away from the sam- 
ple. The dc magnetic field was applied in the direction 
of the incoming neutron beam using a horizontal super- 
conducting magnet. A coil was wound on the sample to 
allow in situ ac magnetic susceptibility measurements. 

Fig.l(a) shows the SANS data at H= 3.0 kOe. The 
Gaussian width data are obtained from fitting the mea- 
sured (Bragg) intensity vs. azimuthal angle to six Gaus- 
sian peaks evenly spaced 60° apart. It is clear that the 
azimuthal widths, a measure of orientational disorder in 
the vortex array, are strongly history dependent. Super- 
cooling and superheating effects are observed for field- 
cooling (FC) and field-cooled-warming (FCW) paths, re- 
spectively. As reported previously |8j, the disordered 
phase at T > T p and the ordered phase at T < T p are 
of their respective thermodynamic ground states. The 
abrupt change in the structure factor 5(g) at the peak 
effect T p depicts a symmetry-breaking phase transition 
from a vortex matter with short-range order to a Bragg 
glass with quasi- long range order p|. The phase tran- 
sition is first order as evidenced by the strong thermal 
hysteresis in S(q). Compared to that at higher fields, 
the metastability region for H = 3.0 kOe is smaller but 
still pronounced. 

We found that the thermal hysteresis of 5(g) observed 
in SANS is strongly field dependent, and the metasta- 
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bility region disappears completely at low fields. Fig.l 
(b) shows the azimuthal width data for H = 2.0 kOc. 
For comparison, the real part x'(T) °f the ac magnetic 
susceptibility is also shown in Fig. 1(b). The dip in 
x'(T) i s a well-established signature of the peak effect 
[l2t Il3l . The history dependence of Bragg-peak width 
is detectable only within 0.1 K of the peak-effect temper- 
ature T p . 

A similar trend is observable in the history dependence 
of the radial widths of the Bragg peaks, as shown in Fig. 2, 
obtained by fitting a single Gaussian function to the in- 
dependence of the SANS intensity. At 3.0 kOe, there 
is a pronounced thermal hysteresis in the radial widths. 
At 2.0 kOe, however, the hysteresis is barely discernable. 
At an even lower field of 1.0 kOe (data not shown), the 
thermal hysteresis in S(q) is undetectable. At H = 1.0 
kOe, a very sharp peak effect (the onset-to-end width = 
40 mK) is still present. Thus we believe the phase tran- 
sition at 1.0 kOe is still first-order but the metastability 
region is too narrow to be resolved in SANS (the tem- 
perature resolution in SANS « 50 mK). Nevertheless, 
the diminishing hysteresis in the Bragg-glass transition 
in the low-field regime suggests that the phase transition 
is becoming continuous and mean-ficld-like. 

The fact that the transition into the Bragg glass is first 
order at high fields, but continuous (mean-field-like) at 
low fields strongly suggests the existence of a multicritical 
point on the phase boundary bordering the Bragg-glass 
on the H — T phase diagram. We show below that this 
multicritical behavior is directly related to the appear- 
ance and the disappearance of the peak effect. 

Fig. 3 shows a three-dimensional plot of the x'(T) as 
a function of temperature and magnetic field in the field 
range of - 5.12 kOe. At high fields, there is a pro- 
nounced peak effect, a characteristic dip in x'{T). With 
decreasing field, the peak effect becomes narrower and 
smaller. For H < 0.8 kOe, there is only a single kink in 
x'(T) corresponding to the mean-field transition H C 2(T). 
According to the x'(T) data in Fig. 3, there is no reentrant 
peak effect at low fields, in contrast to that in 2H-NbSe2 
{TH and YBCO [l6|. The peak effect simply vanishes 
here. 

At higher temperatures above the peak-effect temper- 
ature Tp(H) (or H p (T), used interchangeably), there is a 
smooth step in X '{T). This step, T c3 {H) (or H c3 (T)), de- 
fined in Fig.l (b), is the onset of surface superconductiv- 
ity. The separation between T p and T c3 grows larger with 
increasing magnetic field. Upon cooling, below T c3 (H) 
and towards T p (H), the screening effect in 47rx'(T) in- 
creases gradually. Nevertheless, a less well-defined char- 
acteristic temperature T c2 (H) can be identified to mark 
the onset of bulk superconductivity [see Fig. 1(b) for def- 
inition]. Note that the notation T C 2(H) is used for H > 
0.8 kOe, while H c2 (T) is used for H < 0.8 kOe. 

The results of Fig. 3 are summarized in a new phase di- 
agram of Bragg glass superconductivity in Nb as shown 
in Fig. 4. The measured ratio of H c3 /H C 2 at low temper- 
atures is about 1.60, slightly smaller than the expected 



value of 1.695 by Saint-James and de Gennes |l7|. This 
is likely due to the nonidcal cylinder surface being not 
exactly parallel to the field . The crossing of H c3 (T) and 
H C 2 (T) lines below T c was observed previously ^3 ^| , 
and has been interpreted as due to a depressed BCS gap 
function near the surface ppf. 

The most striking aspect of Fig. 4 is that all four lines, 
H p (T), T c2 (H), H c2 (T), and H c3 (T), meet at a multi- 
critical point (MCP). To determine the nature of a MCP, 
one needs to know how many of these lines are second- 
order phase transitions. In the theory of Saint- James and 
de Gennes 0, T c3 (H) is a continuous phase transition. 
For H < 0.8 kOe, the linear temperature dependence 
of H C 2 (T) follows the expected behavior of a Ginzburg- 
Landau mean- field transition line [2]]]. This is a line of 
continuous phase transitions from the normal state di- 
rectly into an ordered Abrikosov Bragg-glass phase. For 
H > 0.8 kOe, the peak effect T P (H) traces out a line of 
first-order transitions. Across this line, the thermal hys- 
teresis in the structural factor S(q) of the vortex matter 
can be observed, especially striking at high fields. 

The nature of the T C 2(H) line is less clear. The vortex 
matter is liquidlike structurally in the shaded part of the 
phase diagram. Whether this disordered vortex matter 
is a distict thermodynamic phase from the normal state 
is still being debated (see Q and references therein). If 
T C 2{H) is a true second-order phase transition, e.g., as in 
a vortex glass transition [2^, the MCP in Fig. 4 would 
appear to be a tricritical point. On the other hand, we 
found that the measured slopes of the four phase bound- 
aries near the MCP cannot satisfy the requirements |2^| 
imposed on a tricritical point by thermodynamics, but 
are consistent with those for a bicritical point. This leads 
to an important conclusion that only one of the two lines, 
T c2 {H) or H c3 (T), can be related to the MCP. 

For T C 2(H) to be a second-order line for the bicritical 
point, its slope has to be larger (in magnitude) than that 
of H C 2(T), such that the thermodynamics rule [24[ that 
no phase can occupy more than 180° of the phase space 
around a MCP is satisfied. However, this is not observed 
in our data as shown in Fig. 4. For the H c3 (T) line to 
be the relevant one, the ratio of specific-heat jump at 
H c2 (T) over that at H c3 {T) should be 43.6. While this 
large ratio is consistent with the existing specific-heat 
data on Nb [2!|, presently there are no reliable specific- 
heat data near a crossing point of H C 2 (T) and H c3 (T) to 
allow us to make a quantitative comparison. 

We should point out that a similar critical point has 
also been observed in platelet geometries such as MgB2 
[26ll27l ] crystals for which H c3 (T) is not expected to play 
a role in the critical point. In high-T c YBCO, a disap- 
pearance of the first-order transition was also observed 
in the low-field regime, and was interpreted as a critical 
end point |2^. If our critical point in Fig. 4 is also in- 
terpreted as a critical end point, T c2 (i?) would not be a 
true phase transition, and the meeting of H c3 (T) at the 
MCP would be purely coincidental. 

In summary, we found that, in a Nb crystal in which a 
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peak effect in ac magnetic susceptibility and a first-order 
melting (disordering) transition in SANS were found to 
coincide previously, both effects disappear at a low field. 
It is suggested that the appearance or absence of a peak 
effect in a type-II superconductor may be directly corre- 
lated with a multicritical point (MCP) on the Bragg glass 
phase boundary. The existence of a MCP, at which the 
peak effect vanishes, suggests that the origin of the peak 
effect is related to a second-order phase transition at a 
higher temperature. In the sample studied, it appears 
that the MCP may be related to surface superconductiv- 
ity. 
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Figure Captions: 

Figure 1: (a) Temperature and history dependence of 
azimuthal widths of the (1,-1) diffraction peak at Hdc = 
3.0 kOe. The widths are obtained by Gaussian fittings. 
The dashed line is the peak of the peak effect T p at 



this magnetic field based on ac magnetic susceptibility 
measurements. Inset: Experimental configuration, (b) 
Temperature and history dependence of the azimuthal 
widths of (1,-1) diffraction peak at Hdc — 2.0 kOe. 
The ac susceptibility data are also shown for reference. 
Definitions of T P (H), T c2 {H), and T c3 (H) (see below) 
are shown. 

Figure 2: Temperature and history dependence of the 
radial widths of the diffraction peaks at (a) Hd c = 3.0 
kOe, and (b) H dc = 2.0 kOe. 

Figure 3: Three-dimensional (3D) magnetic field and 
temperature dependence of the real part of the ac sus- 
ceptibility 4ttx'(T). Hdc\\H ac . Note that two values of 
ac fields were used in the measurements. For Hd c <3.0 
kOe, H ac = 1.7 Oe and / = 1.0 kHz and for H dc > 3.0 
kOe, H ac = 7.0 Oe and / = 1.0 kHz. The solid and 
dashed lines are guides to eyes. For the ac fields used, T p 
is independent of the ac field amplitude 01 • 

Figure 4: The phase diagram of a weak-pinning Nb 
crystal for the Hd c \\ < 111 > crystallographic direction. 
The upper solid circles (the crosses correspond to mea- 
surements using H ac = 1.7 Oe, see Fig. 3 captions) are 
the peak of the peak effect and the first-order transition 
line; the lower ones are the mean-field transition. The 
open diamonds (two sets are for two values of ac fields, 
see Fig. 3) are H c3 . The multicritical point is indicated 
by the large filled circle. The H c i data (triangle) are es- 
timated from the first penetration in the ac susceptibility 
data. All lines are hand drawn as guides for eyes. 

^ Author to whom correspondence should be addressed. 
Electronic address: xsling@brown.edu. 
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